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ABSTRACT 

Context. For the spectral analysis of high-resolution and high-signal-to-noise (S/N) spectra of hot stars, state-of-the-art non-local 
thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic 
data that is used for their calculation. In a recent analysis of the ultraviolet (UV) spectrum of the DA-type white dwarf G191-B2B, 
21 Zn IV lines were newly identified. Because of the lack of Zn IV data, transition probabilities of the isoelectronic Ge VI were adapted 
for a first, coarse determination of the photospheric Zn abundance. 

Aims. Reliable Zn IV and Zn V oscillator strengths are used to improve the Zn abundance detennination and to identify more Zn lines 
in the spectra of G191-B2B and the DO-type white dwarf RE 0503-289. 

Methods. We performed new calculations of Zn IV and Zn V oscillator strengths to consider their radiative and collisional bound- 
bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of the ZnlV — V spectrum exhibited in high- 
resolution and high-S/N UV observations of G191-B2B and RE 0503-289. 

Results. In the UV spectrum of G191-B2B, we identify 31 ZnlV and 16 ZnV lines. Most of these are identified for the first time 
in any star. We can reproduce well almost all of them at logZn = -5.52 ± 0.2 (mass fraction, about 1.7 times solar). In particular, 
the Zn IV / Zn V ionization equilibrium, which is a very sensitive T e ff indicator, is well reproduced with the previously determined 
T e g = 60000 ± 2000 K and log g = 7.60 ± 0.05. In the spectrum of RE 0503-289, we identified 128 ZnV lines for the first time and 
determined logZn = -3.57 + 0.2 (155 times solar). 

Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed 
ZnlV and ZnV line profiles in two white dwarf (G191-B2B and RE0503-289) ultraviolet spectra were well reproduced with our 
newly calculated oscillator strengths. This allowed us to determine the photospheric Zn abundance of these two stars precisely. 

Key words, atomic data- line: identification - stars: abundances - stars: individual: G191-B2B - virtual observatory tools - 
stars: individual: RE 0503-289 


1. Introduction 

In a recent spectral analysis of the hydrogen-rich DA-type white 
dwarf G191-B2B, Rauch et al. (2013) identified and reproduced 
stellar lines of C, N, O, Al, Si, O, P, S, Fe, Ni, Ge, and Sn. In ad- 
dition, they identified 21 ZnlV lines. The determined Zn abun- 
dance (logarithmic mass fraction of -4.89, 7.5 x solar) was 
uncertain because the unknown Zn IV oscillator strengths were 


* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26666. 

** Based on observations made with the NASA-CNES-CSA Far 
Ultraviolet Spectroscopic Explorer. 

*** Tables 1 and 2 are only available at the CDS via anonymous ftp to 
cdsarc .u-strasbg . fr (13® . 79 . 128 . 5) or via 
http : //cdsarc .u-strasbg . fr/viz-bin/qcat?J/A+A/564/A41 


approximated by values of the isoelectronic Gevi taken from 
Rauch et al. (2012). 

In this paper, we introduce new oscillator strengths for Zn IV 
and ZnV (Sect. 2). Then, we describe briefly our observations 
(Sect. 3), our analysis strategy (Sect. 4), and revisit G191-B2B 
to perform a precise determination of its Zn abundance (Sect. 5). 
The white dwarf RE 0503-289 is hotter than G191-B2B and 
its trans-iron element abundances are strongly oversolar (Werner 
et al. 2012; Rauch et al. 2013) and, thus, it appears promising to 
identify Zn lines. In Sect. 6, we describe our search for these 
and the determination of its Zn abundance. We summarize our 
results and conclude in Sect. 7. 

2. Transition probabilities in Zn IV and Zn V 

Radiative decay rates (oscillator strengths and transition prob- 
abilities) have been computed using the pseudo-relativistic 
Hartree-Fock (HFR) method as described by Cowan (1981). 
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Fig. 1. Grotrian diagrams of our ZnlV (top) and ZnV (bottom) model 
ions. Horizontal bars indicate levels, gray lines represent radiative tran- 
sitions with known / values, respectively. The dashed lines show the 
ionization energies. 


For Zn IV, configuration interaction has been considered among 
the configurations 3d 9 , 3d 8 4s, 3d 8 5s, 3d 8 4d, 3d 8 5d, 3d 7 4s 1 2 * * , 
3d 7 4p 2 , 3d 7 4d 2 , 3d 7 4f 2 , 3d 7 4s5s, 3d 7 4s4d, and 3d 7 4s5d for the 
even parity and 3d 8 4p, 3d 8 5p, 3d 8 4f, 3d 8 5f, 3d 7 4s4p, 3d 7 4s5p, 
3d 7 4s4f, 3d 7 4s5f, and 3d 7 4p4d for the odd parity. Using exper- 
imental energy levels published by Sugar & Musgrove (1995), 
the average energies (E m ), the Slater integrals (F k , G k ), the 
spin-orbit parameters (£ni), and the effective interaction param- 
eters (a, P) corresponding to 3d 9 , 3d s 4s, 3d s 4p configurations 
were optimized using a well-established least-squares fitting 
process minimizing the differences between calculated and ex- 
perimental energy levels within both configurations. In the case 
of Zn V, the configurations included in the HFR model were 3d 8 , 
3d 7 4s, 3d 7 5s, 3d 7 4d, 3d 7 5d, 3d 6 4s 2 , 3d 6 4p 2 , 3d 6 4d 2 , 3d 6 4s5s, 
3d 6 4s4d, 3d 6 4s5d for the even parity and 3d 7 4p, 3d 7 5p, 3d 7 4f, 
3d 7 5f, 3d 6 4s4p, 3d 6 4s5p, 3d 6 4s4f and 3d 6 4p4d for the odd par- 
ity. For this ion, the semi-empirical fitting process was per- 
formed to optimize the radial integrals corresponding to 3d 8 , 
3d 7 4s, and 3d 7 4p configurations using the experimental energy 
levels compiled by Sugar & Musgrove (1995). The HFR oscil- 
lator strengths (log gf) and transition probabilities ( gA , in s -1 ) 
for Zn IV and Zn V spectral lines are reported in Tables 1 and 2, 
respectively, alongside with the numerical values (in cm -1 ) of 
lower and upper energy levels and the corresponding wave- 
lengths (in A). In the last column of each table, we also give the 
cancellation factor CF as defined by Cowan (1981). We note that 
very small values of this factor (typically <0.05) indicate strong 
cancellation effects in the calculation of line strengths. In these 
cases, the corresponding gf and gA values could be very inac- 
curate and so need to be considered with some care. However, 
very few transitions appearing in Tables 1 and 2 are affected by 
these effects. Figure 1 shows Grotrian diagrams of ZnlV and 
Zn V including all levels and transitions from Tables 1 and 2. 



log (m / g/cm 2 ) 

Fig. 2. Ionization fractions of Zn II - VI in our G191-B2B model atmo- 
sphere. m is the column mass, measured from the outer boundary of our 
model atmosphere. The formation depths (i.e., r = 1) of the ZnlV - V 
line cores are marked. 

Table 3. Statistics of our N, O, and Zn model atoms for G191-B2B. 


Ion 

NLTE levels 

LTE levels 

Lines 

N II 

1 

246 

0 

III 

9 

57 

10 

IV 

9 

85 

10 

V 

10 

52 

20 

VI 

1 

0 

0 

O II 

1 

46 

0 

III 

9 

63 

6 

IV 

9 

85 

11 

V 

9 

117 

12 

VI 

10 

75 

19 

VII 

1 

0 

0 

Zn II 

6 

0 

8 

III 

13 

0 

17 

IV 

63 

13 

399 

V 

157 

0 

1878 

VI 

1 

0 

0 


3. Observations 

In this analysis, we use the FUSE 1 spectrum of RE 0503-289 
and the FUSE and HST/STIS 2 spectra G191-B2B that are de- 
scribed in detail by Werner et al. (2012) and Rauch et al. (2013), 
respectively. 

Both FUSE spectra are co-added from all available obser- 
vations of RE 0503-289 and G19 1 -B2B. They cover the wave- 
length range 910 A < A < 1188 A. Their resolving power is 
R = A/AA ~ 20000. The HST/STIS spectrum of G191-B2B is 
co-added from 105 observation with the highest resolution (grat- 
ing E140H, R ~ 1 18 000, 1145 A < A < 1700 A) available via 
MAST. 


4. Model atmospheres and atomic data 

To determine the Zn abundance of G191-B2B, it would be 
straightforward to use the final model of Rauch et al. (2013) 
as well as their model atoms with the only exception that the 
Zn IV and Zn V model ions were replaced by the extended ver- 
sions that consider the newly calculated transition probabilities 

1 Far Ultraviolet Spectroscopic Explorer. 

2 Hubble Space Telescope/Space Telescope Imaging Spectrograph, for 

our high-resolution spectrum of G191-B2B, see http://www. stsci . 

edu/hst/observatory/crds/cal spec . html 
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Fig. 3. Zn IV lines (left panel , marked with their wavelengths in A, blue in the online version) and Zn V lines ( right panel , marked in green) in 
the FUSE (for lines at d < 1 150 A) and HST/STIS (d > 1 150 A) observations of G191-B2B compared with our theoretical line profiles. For the 
identification of other lines, see Rauch et al. (2013). The vertical bar shows 10% of the continuum flux. 
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Table 4. Identified Zn lines in the UV spectrum of G191-B2B. 


Wavelength/A 

Ion Comments 


Theoretical Observed 


1272. 

202 



1272. 

990 

1272.975 

prev. unid. 

1277. 

080 

1277.130 

A A > 0.02 A 

1283. 

478 

1283.525 

Ad > 0.02 A, prev. unid. 

1284. 

721 

1284.740 


1291. 

826 

1291.810 


1296. 

652 

1296.620 


1320. 

704 

1320.725 


1321. 

215 



1322. 

316 

1322.320 


1326. 

774 

1326.735 

Ad > 0.02 A 

1329. 

110 



1330. 

302 

1330.325 


1333. 

180 



1340. 

156 

1340.190 

A/l > 0.02 A 

1342. 

716 

1342.755 

Ad > 0.02 A 

1343. 

750 

1343.815 

Ad > 0.02 A 

1344. 

122 

1344.090 

A/l > 0.02 A 

1347. 

954 

1347.970 


1349. 

876 

1349.895 


1352. 

883 

1352.905 

A/l > 0.02 A 

1356. 

171 

1356.090 

Ad > 0.02 A, blend with 




Zn IV /l 1356.195 A 

1357. 

801 

1357.810 


1359. 

477 

1359.490 


1363. 

432 

1363.420 


1363. 

912 

1363.940 

A/l > 0.02 A 

1365. 

253 

1365.260 


1369. 

510 

1369.515 


1375. 

325 



1377. 

615 

1377.635 


1387. 

694 

1387.720 


1116. 

842 

1116.860 


1120. 

325 

1120.330 


1133. 

031 

1133.060 

Ad > 0.02 A, prev. unid. 

1147. 

020 

1147.040 


1152. 

985 

1152.980 


1155. 

027 

1155.045 


1158. 

759 

1158.750 

prev. unid. 

1174. 

346 

1174.325 

prev. unid. 

1176. 

122 


weak, prev. unid. 

1179. 

969 

1180.005 

A/l > 0.02 A 

1185. 

898 

1185.905 


1185. 

948 

1185.955 


1186. 

057 



1186. 

447 


weak, prev. unid. 

1195. 

745 


weak 

1202. 

906 




Notes. Observed wavelengths are given only in case that they deviate 
from the theoretical wavelengths (cf. Tables 1 and 2). “prev. unid.” de- 
notes lines that were previously listed as unidentified by Rauch et al. 
(2013). 


(Sect. 2). Unfortunately, the employed Tubingen non-local ther- 
modynamic equilibrium (NLTE) model-atmosphere package 
(Werner et al. 2003; Rauch & Deetjen 2003, TMAP 3 ), which 

3 http : //astro . uni - tuebingen . de/~TMAP 



Fig. 4. Theoretical line profiles of the strongest ZnlV lines (left) and 
Zn V lines (right) (marked with their wavelengths from Tables 1 and 2) 
calculated from our model of G191-B2B with a Zn abundance of 
3.0 x 10 -6 (mass fraction) located in the STIS wavelength range com- 
pared with the observation. The lines are shifted to the observation, see 
Table 4. 



Fig. 5. Theoretical line profiles of the strongest ZnV lines calculated 
from our model of RE 0503-289 with a Zn abundance of 2.7 x 10 -4 
(mass fraction) located in the FUSE wavelength range compared with 
the observation. The lines are shifted to the observation, see Table 5. 


is used to calculate plane-parallel, chemically homogeneous, 
metal-line blanketed NLTE model atmospheres, overcharged our 
FORTRAN compilers. The program would not compile if the ar- 
ray sizes were increased further according to the much higher 
number of atomic levels treated in NLTE and the respective 
higher number of radiative and collisional transitions. 

Thus, we decided to reduce the number of N and O lev- 
els treated in NLTE (Table 3) to create a TMAP executable. 
Test calculations have shown that the deviations in temperature 
and density structure between the final model of Rauch et al. 
(2013) and a model with reduced N and O model atoms are neg- 
ligible. Then, the Zn occupation numbers are determined in a 
line-formation calculation, i.e., at fixed temperature and density 
structure. Since Zn opacities were already considered in our start 
model, the atmospheric structure and the background opacities 
are well modeled. 
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Fig. 6. ZnV lines in the FUSE observation of RE 0503-289 compared with our theoretical line profiles. 


All model atoms (including Zn) are provided via the 
Tubingen Model-Atom Database (TMAD 4 , Rauch & Deetjen 
2003), that has been set up within a project of the German 
Astrophysical Virtual Observatory (GAVO 5 ). All SEDs that 
were calculated for this analysis are available via the registered 
Theoretical Stellar Spectra Access (TheoSSA 6 ) VO service. 

4 http : //astro . uni - tuebingen . de/~TMAD 

5 http://www.g-vo.org 

6 http://dc.g-vo.org/theossa 


5. The photospheric Zn abundance in G191-B2B 

Zn IV and Zn V are the dominant ionization stages of Zn in the 
atmosphere of G191-B2B (Fig. 2). Therefore, we closely in- 
spected the available spectra for lines of these ions. 

In the FUSE and HST/STIS observations of G191-B2B (cf. 
Rauch et al. 2013)we identified31 Zniv (lOnew identifications) 
and 16 ZnV (all new) lines. The observed wavelength positions 
(a radial velocity of u ra( j = 22.1kms _1 was applied according 
to Holberg et al. 1994; Rauch et al. 2013) deviate partly from 
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Table 5. Like Table 4, for RE 0503-289. 


Table 5. continued. 


Wavelength/A 

Comments 

Ion 

Wavelength/A 

■ Comments 

Theoretical 

Observed 

Theoretical 

Observed 

Zn v 1029. 992 




1138. 

248 



1035. 859 




1138. 

497 



1035. 887 




1139. 

278 

1139.220 

Ad > 0.02 A 

1043. 353 







blend with Ge V 

1055. 878 




1140. 

703 



1063. 979 




1141. 

003 

1141.015 


1068. 284 




1141. 

095 



1069. 674 


blend with 0 V, Ga V, Ge V 


1141. 

344 



1069. 764 




1142. 

792 



1071. 501 




1142. 

925 



1072. 992 

1072.950 

A A > 0.02 A 


1142. 

938 


prev. unid. 

1074. 241 

1074.265 

A d > 0.02 A 


1143. 

196 



1075. 171 

1075.050 

Ad > 0.02 A 


1143. 

403 



1076. 239 




1144. 

136 

1144.160 

Ad > 0.02 A 

1076. 878 

1076.895 



1145. 

151 



1090. 831 

1090.800 

Ad > 0.02 A 


1146. 

057 



1094. 088 

1094.110 

blend with N IV 


1146. 

149 


too strong in model 

1095. 774 


blend with Ge V 


1147. 

020 

1147.040 


1095. 797 




1147. 

371 

1147.425 

Ad > 0.02 A 

1095. 961 

1095.945 



1148. 

922 

1148.915 


1109. 078 

1109.110 

Ad > 0.02 A 


1149. 

398 

1149.370 

Ad > 0.02 A 



blend with C IV 










1149. 

486 



1109. 166 


blend with C IV 






1110. 821 

1110.810 

weak 










1149. 

873 

1149.855 


1111. 530 


blend with C III, O IV 










1 150. 

743 















1151. 

368 



1112. 829 








1114. 482 




1152. 

985 

1152.980 


1115. 266 

1115.295 

Ad > 0.02 A 


1 153. 

160 



1115. 668 

1115.695 

Ad > 0.02 A 


1155. 

027 

1155.045 


1116. 630 


too strong in model 


1155. 

725 

1158.750 


1116. 842 

1116.860 



1156. 

520 



1119. 950 

1119.940 



1156. 

885 



1120. 101 

1120.080 

Ad > 0.02 A 


1158. 

122 



1120. 325 




1158. 

475 



1121. 109 

1121.095 



1158. 

759 

1158.750 


1121. 524 


weak 


1160. 

091 


weak 

1122. 502 


blend with Si IV 


1160. 

221 



1123. 127 


blend with Ga V 


1161. 

971 



1125. 019 

1125.050 

Ad > 0.02 A 


1162. 

281 



1125. 048 

1125.060 



1162. 

401 



1129. 898 


blend with Ga V 


1164. 

090 


blend of 

1130. 051 







Znv dd 1165.082, 1164.101 A 

1130. 242 




1164. 

632 


blend with O IV 

1131. 242 

1131.250 

prev. unid. 


1165. 

189 



1131. 788 




1165. 

706 


blend with C III 

1131. 863 




1165. 

716 


blend with C III 

1132. 271 

1132.290 



1165. 

880 


blend with Xe VI 

1132. 659 


blend with N IV 


1170. 

105 



1133. 031 

1133.060 

Ad > 0.02 A 


1171. 

106 

1171.130 


1133. 128 




1173. 

366 



1133. 278 

1133.300 

Ad > 0.02 A 


1174. 

346 

1174.325 


1133. 498 




1174. 

945 


blend with C III 

1135. 324 




1176. 

122 


blend with C III 

1135. 588 




1178. 

759 



1136. 603 




1179. 

145 



1136. 986 

1137.000 



1179. 

179 



1137. 625 




1179. 

969 

1180.005 

Ad > 0.02 A 
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Table 5. continued. 


Wavelength/A 

Ion Comments 

Theoretical Observed 


1180. 018 


blend of 

Znv AA 1 180.018, 1 180.025 A 

1182. 019 

1182. 567 

1183. 041 
1183. 100 
1183. 158 
1183. 314 
1185. 619 
1185. 645 
1185. 676 
1185. 898 

1185.905 


1185. 948 

1185.955 

blend of 

Znv AA 1 185.948, 1 185.961 A 

1186. 057 
1186. 447 

1186.420 

Ad > 0.02 A 

1187. 664 


too strong in model 

1187. 706 


too strong in model 

those given in Tables 1 and 2 by 

some hundredths of an A. The 

good agreement of the strongest, unshifted lines in our model 7 


(Fig. 3) with the observations permits to shift the lines to ob- 
served absorption features in their closest vicinity. The reason 
for this uncertainty is most likely the limited accuracy of the 
Zn IV and Zn V energy levels from which the wavelengths of the 
line transitions were calculated. The identified lines are summa- 
rized in Table 4. 

Our calculations have shown that the Zniv/Zn V ionization 
equilibrium at r e ff = 60000K and log g = 7.6 (cf. Rauch et al. 
2013) is well reproduced (Figs. 3, 4). On the other hand, the 
Zn abundance given by Rauch et al. (2013, 1.3 x 10 -5 ± 0.5 dex 
by mass) is too high. We reduced it to 3.0 x 10 -6 ± 0.2 dex 
(about 1.7 times solar, following Asplund et al. 2009) to repro- 
duce the observed Zn lines best. This agrees with the previous 
value within the error limits. Even a solar Zn abundance, how- 
ever, is possible within the error limits. 

6. The photospheric Zn abundance in RE 0503-289 

Our inspection of all UV spectra that were already used by 
Werner et al. (2012) has shown that only RE 0503-289 exhibits 
prominent Zn lines (Fig. 5). We find that in its FUSE spectrum 
(»md = 26.3kms _1 ), a rich Znv spectrum of 128 lines (Fig. 6, 
Table 5) is present. The synthetic spectrum of our final model 
shows more, weak lines that do not have an unambiguous line 
identification due to the signal-to-noise ratio (S/N) of the ob- 
servation. The model’s prediction of their relative line strengths 
facilitates to distinguish between noise and “real” lines in the 
observation and, hence, to identify even such weak lines. In gen- 
eral, all lines with oscillator strengths gf > 0.01 can be de- 
tected. 

Dreizler & Werner (1996) determined T e s = 70 000 ± 4000 K 
and log g = 7.50 ± 0.25 for RE 0503-289. This was recently 
verified by well-matched ionization equilibria of Kr and Xe 

7 All synthetic spectra shown in this paper are convolved with 
Gaussians to match the spectral resolution (FUSE: FWHM = 0.06 A, 
STIS: FWHM = 0.01 A). 



Fig. 7. Like Fig. 2, for our RE 0503-289 model atmosphere. 

(Werner et al. 20 1 2, Kr VI / Kr VII, Xe VI / Xe VII) and Ge (Rauch 
et al. 2012, Gev/Gevi). We adopt these values and start our 
calculation based on the final model of Werner et al. (2012) that 
considers opacities of He, C, N, O, Ge, Kr, Xe, and of the iron- 
group elements (Ca - Ni). We follow the same strategy described 
in Sect. 4, but this time, we reduced the size of the Ge model 
atom in the line-formation calculations. 

The higher T t g compared to that of G191-B2B shifts the 
Zn ionization equilibrium strongly towards higher ionization 
(Fig. 7). Znv remains dominant while Zniv is less occupied 
by a factor of about 100 at all depths. The ionization fraction of 
Zn VI is also much below that of Zn V and we only expect weak 
lines. The strongest Znvi lines 8 9 are located in the soft X-ray 
to EUV g wavelength range where we do not have high-quality 
observations to evaluate. 

We determine a Zn abundance of 2.7 x 10 -4 ± 0.2 dex 
(about 155 x solar) to reproduce the observed Znv line profiles 
best (Fig. 5). 

7. Results and conclusions 

The identified Zniv and Znv lines in the high-resolution 
UV spectra of G191-B2B and RE 0503-289 are well repro- 
duced with our newly calculated oscillator strengths by our 
NLTE model-atmosphere calculations. 

We determined photospheric abundances of log Zn = 
-5.52 ± 0.2 (mass fraction, 1.9-4. 8 x 10 -6 , 1. 1-2.8 times the 
solar abundance) and log Zn = -3.57 ± 0.2 (1.7-4. 3 x 10 -4 , 
98-248 times solar) for the DA-type white dwarf G191-B2B 
and the DO-type white dwarf RE 0503-289, respectively. The 
highly supersolar Zn abundance is in line with the high abun- 
dances of trans-iron elements Ge (650 x solar, Rauch et al. 
2012), Kr (450 x solar), Xe (3800 x solar, Werner et al. 2012) 
in RE 0503-289. 

The identification of new lines due to trans-iron elements, 
e.g., Ga, Ge, As, Se, Kr, Mo, Sn, Te, I, and Xe (Werner 
et al. 2012) and Zn (Rauch et al. 2013, and in this paper) in 
G191-B2B and RE 0503-289 promises to help enhance the un- 
derstanding of extremely metal-rich white dwarf photospheres 
and their relation to AGB and post-AGB stellar evolution. Their 
reproduction, i.e., the precise abundance determination, e.g., 
of Kr and Xe (Werner et al. 2012), Ge and Sn (Rauch et al. 
2012), and Zn (this paper) is strongly dependent on the available 
atomic data. This remains a challenge for atomic and theoretical 
physicists. 
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